Metacinnabar (/3-HgS): a strong 3D topological insulator with highly anisotropic 

surface states 
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We establish the presence of topologically protected edge states on the (001) surface of HgS in the 
zinc-blende structure using density-functional electronic structure calculations. The Dirac point of 
the edge state cone is very close to the bulk valence band maximum. The Dirac cone is extremely 
anisotropic with a very large electron velocity along one diagonal of the surface elementary cell x' 
and a nearly flat dispersion in the perpendicular direction y' . The strong anisotropy originates from 
a broken fourfold rotoinversion symmetry at the surface. 



PACS numbers: 

Motivated by the ongoing search for spintronics mate- 
rials, recent theoretical work on the effect of spin-orbit 
interactions on the band structure of solids predicted the 
existence of novel topological insulators (TPs) in two [T|- 
[3] and three dimensions [1H7] . Topological insulators dis- 
tinguish themselves from ordinary insulators due to the 
presence of a non-trivial topological invariant in the bulk 
band structure [IHE]. Although there is an excitation 
gap in the bulk as in conventional band insulators, the 
topological bulk properties dictate that there be gapless 
modes at the interface between a TI and a topologically 
trivial state, for instance the vacuum. The surface states 
of a three-dimensional (3D) TI form massless Dirac cones, 
with a single Dirac cone being the simplest case. The 
potential applicability for spintronics arises because the 
surface states of such a 3D TI appear in time-reversed 
pairs in which the two electrons have both opposite spin 
and velocity. 

From the material's viewpoint, two-dimensional (2D) 
TI behavior was experimentally observed in HgTe quan- 
tum wells [3] , soon after the theoretical prediction [3] . 3D 
TPs discovered so far are the binary compounds Bi 2 Se 3 , 
Sb 2 Te 3 , Sb 2 Te 3 QME1] and the alloy Bii^Sb^ g] H3] 
and very recently, 3D-TI behavior was proposed in cer- 
tain ternary Heusler compounds [TU [15]. In the for- 
mer the presence of Dirac cones has been established by 
photo-emission, but direct electronic transport evidence 
for metallic surface states on top of an insulating bulk 
is lacking, mainly because of the small bulk gaps. It is 
interesting to note that before the advent of TI's, in semi- 
metallic bismuth surface states were identified that are 
in hindsight very similar to the chiral surface states of 
TI's [TS] . 

Here we establish that metacinnabar, the binary com- 
pound /3-HgS, is a strong 3D topological insulator. As 
opposed to other TI's discovered to date, the topologi- 
cally protected Dirac cone on the (001) surface of HgS in 
the zinc-blende structure is extremely anisotropic. This 
originates from a broken fourfold rotoinversion symmetry 



at the surface. The presence of strongly anisotropic topo- 
logically protected modes implies highly one-directional 
spin and electron transport properties providing a natu- 
ral route to create potentially ideal quantum wires. As 
this electronic anisotropy is absent in the bulk because it 
is fourfold rotoinversion invariant, the presence of signif- 
icant spin and charge transport anisotropy at the (001) 
surface of HgS will be a smoking gun for its TI nature. 

Metacinnabar, /3-HgS, is a gray-black natural mineral, 
found in mercury deposits that are formed near-surface, 
under low-temperature conditions. The bulk electronic 
structure of HgS in the zinc-blende structure has been 
considered before, with contradicting conclusion on the 
existence of a band gap inversion [T7HI9] . Delin [17] and 
Cardona et al [18] find HgS to be insulating with an un- 
conventional band order at the T point: the conduction 
band minimum has p-character and the s-band, which 
usually builds the conduction band minimum, is about 
one eV below the Fermi level. However, according to the 
electronic structure calculation of Moon and Wei |19) . 
such a band inversion due to spin-orbit interaction is ab- 
sent. For lattice structures with inversion symmetry the 
parity criteria of Fu et al [5] guarantee that an inverted 
semiconductor is a topological insulator, but owing to 
the lack of inversion symmetry of the zinc-blende crystal 
structure of /3-HgS these criteria do not apply. We will 
therefore establish the TI nature of /3-HgS via the bulk- 
boundary correspondence and explicitly calculate its sur- 
face states. 

To this end we performed high precision, all electron, 
full potential and fully relativistic electronic structure 
calculations for /3-HgS using the FPLO (full potential 
local orbital) code, version 9.01-35 [3U][ST]. We unam- 
biguously establish the presence of an inverted band gap 
in the bulk material. Investigating in detail the (001) 
Hg-terminated surface in a slab geometry we demon- 
strate the presence of a highly anisotropic Dirac cone. 
This proves that HgS is a topological insulator. The 
anisotropy of the cone implies that the velocity and mo- 
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FIG. 1: Band structure of HgS in the bulk. It is an inverted 
semiconductor which is visible at the F-point: small gap and 
electron-like dispersion at the top of the valence band. 

bility of the charge carriers at the surface is very high in 
one direction and low in the direction perpendicular to 
it. Thus it will strongly affect surface spin and charge 
transport. 

To calculate the bulk band structure of the zinc-blende 
crystal we used the experimental lattice parameter ao = 
5.85 A, incidentally equal to the theoretical value (Fig.[l|. 
Without spin-orbit (SO) coupling, all three 3p-orbitals of 
S would be degenerate leading to metallic behavior. The 
SO-coupling splits the 3p-manifold into a lower quartet 
(3p 3 / 2 ) and an upper doublet (3p!/ 2 ) with a gap of 109 
meV at the T-point. The indirect gap between conduc- 
tion and valence band is 42 meV. The band of s-character 
(here Hg 6s) which builds the conduction band minimum 
in ordinary semiconductors is found at 0.65 eV below the 
Fermi level. This level inversion, in connection with the 
electron like dispersion of the valence band at the T-point 
implies that HgS is an inverted semiconductor, which is 
a prerequisite for it to be a topological insulator. 

The insulating character of HgS is related to the fact 
that the sulfur quartet (3^3/2) is below the doublet 
(3^1/2) which contradicts the third Hund's rule for pure 
p-orbitals. The related compounds HgSe and HgTe, in 
which the Hg 6s level is far below the Fermi level too, 
have the "correct order" with the doublet below the quar- 
tet. Accordingly, there is no gap in those two compounds. 
The reversed SO-order in HgS (quartet below doublet) 
comes about because of a small but significant contri- 
bution of Hg 5d orbitals whose SO-coupling dominates 
the one of S 3p and reverses the sign. So, contrary to 
naive expectations it is actually the smallness of the SO- 
coupling parameter of the sulfur p orbitals (in contrast 
to the p orbitals of Se or Te) [17] which allows HgS to 
become a topological insulator. 

Since HgS has an inverted level order at the T-point, its 
topological invariant can be non-trivial and thus be differ- 
ent from the vacuum. This implies a change of topology 
at any surface and hence the presence of gapless surface 
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FIG. 2: Schematic representation of the (001) surface of HgS. 
Shown are the mercury top layer (full circles: Hg atoms, 
dashed circle: surplus Hg atom that has been removed) and 
the sub-surface sulfur layer (squares: S). At the surface the 
rotoinversion symmetry of the bulk crystal structure is bro- 
ken. 

states. An explicit calculation of the topological invari- 
ant of HgS is technically more involved as the simple 
definition for lattices with inversion symmetry cannot be 
applied for zinc-blende crystals which lack this symme- 
try [22]. Therefore, we studied in detail the (001) surface 
of /3-HgS and calculated its electronic structure explicitly, 
establishing the topological structure of the bulk via the 
bulk-boundary correspondence [7] [8] . To this purpose we 
performed electronic structure calculations of HgS slabs, 
varying the thickness. To guarantee electro-neutrality of 
the surface each second Hg-atom has to be removed from 
a Hg-terminated surface (see Fig. [2]). The surplus Hg- 
atoms in the top and bottom layers were removed in a 
checker-board pattern, so that the slab has tetragonal 
symmetry. Such a surface reconstruction is known from 
similar surfaces like for instance the (100) surface of WO3 
|23j . We have also considered S-terminated surfaces but 
they are less relevant to the present discussion as we find 
such a termination to be 265 meV per surface elemen- 
tary cell higher in energy and therefore substantially less 
stable. 

The HgS slab that we consider has a tetragonal su- 
percell (space group 111) consisting of (001) stacks of 
conventional cubic elementary cells, spaced by an empty 
(vacuum) layer of 23.4 A thickness. The experimental 
cubic lattice parameter, ao = 5.85 A, was used as in 
the bulk calculations. Self-consistent calculations were 
performed in a scalar-relativistic mode, using the local 
density approximation in the parameterization of Ref. 
1241 The presented band structures were evaluated in 
the fully rclativistic mode (four-component Dirac-Kohn- 
Sham scheme) using a dense fc-point mesh of 25 x 25 x 1 
intervals in the full Brillouin zone and linear tetrahedron 
integration with Blochl corrections in order to ensure the 
correct position of the Fermi level. The valence basis set 
comprised sulfur (2s, 2p, 3s, 3p, 3d, 4s, 4p) and mercury 
(5s, 5p, 5d, 6s, 6p, 6d, 7s) states. 

We systematically studied slabs containing n elemen- 
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FIG. 3: (Color) Band structure of the 8-layer slab (lines) 
superimposed with the bulk band structure projected onto 
the surface (shaded area). The relevant surface bands are 
indicated in red (top surface) and blue (bottom surface). 



tary cells with n ranging from 3 to 8. With increas- 
ing slab thickness the gap caused by the coupling of the 
surface states at opposite sides of the stack gradually 
closes (see Fig. [4] a) and a band crossing appears which 
is most clearly seen for the 8-layer slab the band struc- 
ture of which is shown in Fig. [3] In this figure the band 
structure of the slab calculation is superimposed with the 
bulk band structure projected onto the surface. There 
clearly are four bands crossing the Fermi level at the T- 
point. Two of them belong to the top and two of them 
to the bottom layer, i.e., there is one Dirac cone per sur- 
face layer. Therefore HgS is a strong TI. The calculated 
weights of these bands clearly show that they are surface 
bands with predominantly 3p character from the sulfur 
atoms in the sub-surface atomic layer. It is remarkable 
that these surface states decay rather slowly into the bulk 
(see Fig. [4] b) , which explains the remaining (tiny) gap 
for finite layer thickness. 

The two resulting Dirac cones are almost degenerate in 
the direction T-X but in the direction T-M they split up 
because of the broken fourfold rotoinversion symmetry at 
the surface. The direction T-M corresponds to the diag- 
onal k x > = (k x + fcj,)/v / 2 which is identical to the other 
diagonal k y > — (k x — k y )/\/2 in the tetragonal, rotoin- 
version symmetric electronic structure of the slab (space 
group 111), but not at either of the surfaces. In fact, the 
Dirac cone with the large dispersion along k x * belongs to 
the top surface and the other one with two nearly flat 
bands above and below the crossing point corresponds 
to the bottom surface. Along the other diagonal, these 
roles are interchanged. Indeed, considering just the top 
Hg-layer and the sulfur atoms in the sub-surface layer, 
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FIG. 4: (a) Width of the gap as function of the number of 
elementary cells (EC) and (b) sum of 3p weights over the four 
surface states at T, resolved for each sulfur layer in the slab 
of 8 EC. 



one observes a chain-like structure along one of the diag- 
onals (see Fig. |2|. As will become more clear from the 
continuum model below, the anisotropy can effectively 
be described by a crystal field splitting of p-orbitals that 
build the Dirac cone of the surface. 

To illustrate our band structure results we construct 
a simple two-dimensional (2D) model that describes the 
anisotropic "Dirac cone". This model displays essential 
features of the slab band structure in the neighborhood 
of the T-point. It uses two perpendicular p-like orbitals 
with different on-site energies e x > and e y i (crystal field 
splitting) and directed along the two diagonals of the el- 
ementary cell x' — (x + g/)/v2 and y' — (x — y)/>/2. 
The crystal field splitting accounts for the fact that the 
fourfold rotoinversion symmetry is broken at the surface. 
The two p-like orbitals will be denoted as \p x i s) (for the 
two spin directions s =t or \) and \p y > s) . The p-like or- 
bitals are no eigenstates of the local SO-coupling operator 
Asl but the corresponding coupling elements are easily 
calculated. The SO-coupling gives also rise to spin-flip 
hopping in both directions a,k x iOi and uk y i<j\. Finally, 
there are the known tight-binding hoping elements corre- 
sponding to (T-bonds (strength A) and 7r-bonds (strength 
B). In the continuum limit this Hamiltonian is cast in 
the form of the 4 by 4 matrix 



H = 



E x ,<7 + ak x >a 2 -i(X/2)a 3 
i (A/2) (73 E y '(jo + a.k y io\ 



(1) 



where o~o and 01,2,3 are the identity and three Pauli ma- 
trices, respectively, and E x > — e x > + Ak 2 x , + Bk 2 , and 
E v , = e y ,+Ak*,+Bkl,. 

All momenta k x > and k y > are expressed in units of l/oo 
with do — 5.85 A. The two upper bands of that model fit 
the band structure of the 8-layer slab in the neighborhood 
of the T-point with the following parameters (see Fig. [5| : 
A = -0.05, a = 0.25, AE = e x ,-e y > = 0.1, A = 0.8, and B= 
-0.1 (all values in eV). For the two relevant surface bands, 
in the neighborhood of the F-point, the model ([T]) can be 
further reduced to the anisotropic, massless Dirac Hamil- 
tonian H = a X 'k x iai + a y ik y r<T2, with a x > = a cos 2 7, 
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FIG. 5: (Color) Schematic representation of the anisotropic 
Dirac cone (top) emerging from the surface bands of an 8-layer 
slab of HgS (below). The band structure data near the Fermi 
level are shown as points and the fit with the effective long- 
wavelength model as lines. The points in the surface Brillouin 
zone are denoted as X — (ir/ao,0), M — (ir/ao, ir/ao), and 



(n/ao, —n/ao) in the (k x , k y ) coordinate system. 



a y r = a sin 2 7, and tan(27) = \\\/AE. From our fit we 
find the anisotropy of the Fermi velocities to be v x i /v y i 
= 1/ tan 2 7 ps 18, which is very large indeed. An addi- 
tional refinement of this number for the anisotropy ratio 
requires further density-functional calculations on thicker 
slabs. If we neglect in the model description the crystal- 
field splitting, the Dirac cone becomes again isotropic, as 
it is in for instance Bi 2 Se3. 

To summarize, our density-functional calculations 
identify metacinnabar, /3-HgS, to be a topological in- 
sulator. A detailed study of the (001) surface reveals 
that its Dirac cone is very anisotropic. The velocity 
anisotropy v x i/v y i for an eight unit cell slab is ~18. 
This strong anisotropy of the topological electronic sur- 
face states arises from a broken fourfold rotoinversion 
symmetry and sharply contrasts the perfectly isotropic 
bulk band structure. From this observation we expect 
instead an isotropic Dirac cone at the (111) surface since 
the threefold rotation axis along the (111) direction of 
the bulk remains unbroken at the (111) surface. These 
predictions can directly be tested by experimental tech- 
niques that probe surfaces states, for instance by high 
precision photo-emission experiments or STM. In partic- 
ular spin and charge transport at the (001) surface will be 
very sensitive to the surface anisotropy of the electronic 



structure that is absent in bulk. 
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